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SUMMARY 


Eight alternative modeling techniques to specify prescribed temperature at grid or 
scalar points for transient thermal analyses are presented. Four cases are for con- 
stant temperatures, and the others are time-varying temperature functions. Theoreti- 
cal explications and detailed listing of input data cards used for illustrating different 
modelings are given. It is shown that the NTA can be exploited tp extend beyond its 
normal capabilities through innovative modeling techniques. 

In addition, the effect of node valency on the energy distribution to grid points is 
illustrated and discussed. Guidelines to delineate this effect are given. 


INTRODUCTION 


The NASTRAN Thermal Analyzer (NTA) is a finite-element based general purpose 
heat transfer computer program integrated in the NASTRAN system (refs. 1 and 2). Its 
development was intended to provide anal 3 d;ical analysis capabilities in the multiple dis- 
ciplinary areas, fecial attention was directed to the interface problem which inter- 
feres with reliable predictions of the thermo-stresses or deflations, especially for 
very large sized structures and complex configurations. The NTA was developed by 
making an effective use of the functional modules available in NASTRAN to their fullest 
extent. In addition, new modules comprising new elements and solution algorithms to 
satisfy the unique requirements posed by thermal applications were added. The NTA is 
indeed very versatile and flexible. Analysis capabilities, however can be further ex- 
ploited by innovative modeling techniques to compensate for some of the apparent re- 
strictions of the program, which were not implemented solely for economic reasons. 
Working around such restrictions through intelligent modeling is an achievable practice 
in engineering applications. 

This paper deals with two essential subjects relevant to the NTA modeling. The 
problems considered are: 
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(1) Alternative modeling techniques for specifying the prescribed temperatures at 

grid or scalar points in transient thermal analyses 

A total of eight distinct modeling techniques are presented. Four cases concern 
with a constant temperature value, and the others, a time- varying temperature function. 
For the ease of reference and comparison, all cases are illustrated using a single 
basic physical problem. Details of altering the input data cards from case to case are 
given. Theoretical explications which amplify the basis of modeling techniques are in- 
cluded as necessary. Accuracy of solutions along with individual merits and disadvan- 
tages is discussed. While not attempting to be exhaustive, these illustrations will pro- 
vide general NTA users with a practical guide supplementing what has been documented 
and demonstrated in the two-volume set manuals (refs. 3 and 4). 

(2) The effect of the node valency on energy distribution to grid points 

The energy distribution to grid points attributed to the effect of the node valency 
is examplified by a unit surface of a rectangular area which is heated by an external 
heat flux. This rectangular surface is represented in six different models using two- 
dimensional triangular and quadrilateral elements. Users are reminded by facts to be 
cautious in a modeling process. Guidelines and suggestions to delineate this effect are 
given. 


MODELING TECHNIQUES FOR PRESCRIBED TEMPERATURES 


Unlike the steady-state linear and nonlinear cases (using the rigid formats APP 
HEAT, SOL 1 and SOL 3, respectively), for which the NTA permits direct modeling of 
the prescribed temperature at grid or scalar points via the provided single-point- 
constraint ( SPCl *) and temperature (TEMP and TEMPD) cards directly, the modeling 
to specify the prescribed temperatures in both linear and nonlinear transient thermal 
analyses (using the rigid format APP HEAT, SOL 9) requires indirect and more intri- 
cate techniques. A well adopted method to simulate the effect of a prescribed tempera- 
ture Tj (i denotes grid or scalar point identification) is to connect a good thermal con- 
ductor between the point i and a datum of temperature at zero degree. The thermal 
conductor has a very large thermal conductance, where the subscripts i and o 
denote the terminals of connection. A power or rate of heat flow Q. is then applied to 
that point to yield the desired value of a constrained temperature, according to 


♦The names of actual NTA cards are capitalized and underlined. 
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where = 0 has been implied, and is a thermal conductance defined by 


K . 

1 o 



i o 


( 2 ) 


k, A and £ are thermal conductivity, cross-sectional area and length of the conductor, 
respectively. For instance, if a constant value of 512'’C is prescribed at the point i, a 
conductor with a very large thermal conductance, say = 10® W/°C is chosen, and a 
power, Qi = 512 x 10® W, is then applied to i, the desired value of the prescribed tem- 
peratime is obtained from equation (1), i.e. 


512 X 10® 
T. = 


512°C 


( 3 ) 


This method is due directly to an analogy be^een the thermal and structural systems, 
which is depicted in figure 1. A train of the NTA Bulk Data cards consisting of 
CELAS2 (to enter and link between i and o), DAREA (to enter Q^) and TLOAD2 (to 
apply Q j to i and specify the transient loading conditions) is generally employed. 

Based on the same principle, the modeling technique for a constant prescribed 
temperature at a grid point can be extended to the case of time-varying temperatures. 
This can be accomplished by a substitution of a time-dependent thermal load Q. (t) for 
Q. , a constant value, in equation (1), provided admissible quantities and connectable 
NTA Bulk Data cards are selected properly and referenced c^rectly. This approach 
is, by no means, the only modelii^ technique available to the NTA users. Ei^t options 
of modeling techniques to specify constant and time-varying temperatures at grid or 
scalar points will be systematically presented. All cases will be illustrated by a single 
simple physical problem which does not lose any generalities insofar as the thermal 
boundary conditions are concerned. Consequently, any alterations made to the input 
data cards for individual cases can be easily identified and compared, and the accuracy 
of solutions easily evaluated. 

The physical problem is that of determining the transient temperature response 
of a composite slab of infinite extent during a 45-second period (see fig. 2). The rear 
surface of the insulation layer, having the value of e = a = 0.9, dissipates heat radi- 
atively to an environment of 0°C. The front surface of the metallic slab in contact with 
a hot flowing gas has a convective film coefficient of h = 0.071 W/cm^-^C. The thermo- 
physical properties and dimensions are tabulated in table 1. A uniform initial temper- 
ature of 20®C throughout the composite slab is assumed. Temperatures of the flowing 
hot gas are specified as; 


(1) A constant prescribed value maintained at = 1927 °C, 

(2) A time-varying temperature function Tj (t) which is given in figure 3, 

One-dimensional rod elements are sufficient to represent the described problem 
in a finite-element model. A sketch of the composite slab and its finite- element 
representation is shown in figure 2. 

A listing of that portion of the input data deck, which describes the problem in the 
NTA model and is common to all eight cases is shown in figure 4. With comments 
appropriately added, functions of various sections of the Bulk Data Deck are self- 
explanatory, The blank space boimded by the two lines of dotted stars is reserved for 
inserting different packets of the NTA Bulk Data cards which model the prescribed 
temperatures for the hot flowing gas that is represented by GRID or S POINT 200.* 

A description and discussion for each individual case is in order. Relevant Bulk Data 
cards will be listed. 


For a Constant Prescribed Temperature 

(1) To specify a constant temperature via CELAS2 , DARE A , and TLOAD2 

As described previously, this method is a common modeling technique. The Bulk 
Data cards which fill the blank space reserved in the input data deck (fig. 4) are shown 
in figure 5, Case (1) of the attached listing of computer printout excerpts. 

A fictitious large thermal conductance = 10*’ W/cm-°C is entered in the 3rd 
field of CELAS2 . Its 4th field identifies the point which is numbered 200 and is of the 
SPOINT type as indicated by the 0 in the 5th field. The two blanks in the following two 
fields imply that the other end of the thermal conductor is connected to a datum of 0°C. 
The quantity of power, = 1927 x 10® W, is entered in the 5th field in the thermal 
load card DAREA , whose identification 501 appearing in its 2nd field is referenced by 
TLOAD2 in the latter's third field. The set identification 401 of the TLOAD2 must be 
referenced by a DLOAD card in the Case Control Deck to activate the selection of an 
applied thermal load in the transient solution algorithm. Other filled data in these Bulk 
Data cards are obvious when the NTA manual (ref. 3) is consulted. No repetition need 
be made. 

(2) To specify a constant temperature via CELAS2 , SLQAD , and TLOAD2 

To illustrate how a variety of other available thermal load cards may be substi- 
tuted for DAREA in the preceding case, SLQAD is selected for its simplicity. The 


*Unless otherwise noted to use GRID 200 as a restriction, GRID 200 and SPOINT 200 
are interchangeable in all other cases where SPOINT 200 is explicitly indicated. 
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packet of Bulk Data cards for modeling a constant temperature at the S POIN T 200 is 
shown in figure 5, Case (2). Results are identical to that of the preceding case. Other 
thermal load cards such as QVECT , QBDYl, QBDY2 , and QURDY may also serve as 
the replacement for DAREA , but generally require the use of associated cards as an 
inherent requirement. Later cases will illustrate this point. 

(3) To specify a constant temperature via CDAMP2 

The apparent attractiveness of this modeling technique is attributed to its model- 
ing simplicity. A CDAMP2 card loaded with an extremely large thermal mass would 
maintain a prescribed temperature at the applied point, as long as sufficient amount 
of thermal power is supplied. While achieving the intended boundary condition satisfac- 
torily, the addition of a fictitious, very large thermal mass as a component in the cor- 
responding matrix would produce an unstable solution. As a result, accuracy of the 
solution is questionable. Furthermore, there are restrictions associated with this 
technique: Only GRID is permitted to be referenced by CDAMP2 for addii^ a thermal 
mass. A restriction (defect) of the program code requires the presence of CELAS2, 
DAREA , and TLOAD2 in the input data deck even though they would not contribute to the 
result of solution. Although the use of this technique is not recommended. Bulk Data 
cards used in modeling are shown in figure 5, Case (3) for reference. 

(4) To specify a constant temperature via E POINT , DMIG , DAREA and TLOAD2 

The demonstrated modeling techniques in the preceding cases are incapable of 
providing information on the amount of the heat of constraint required to sustain the 
prescribed temperature at a point. In each case, a very large fictRious thermal load 
or a thermal mass was applied to achieve the desired effect. The NTA permits a direct 
input to the conduction matrix through DMIG which can be utilized together with 
E POINT , DAREA and TLOAD2 to enter a true value of the prescribed temperature and 
yet yield the value of the heat of constraint as a part of the solution. Some manipula- 
tions are required, however, to accommodate this scheme. It's theoretical basis will 
be discussed first. 

In transient thermal analysis, the general differential equation in matrix form is 


[C] {f} + [K] {T} = {Q^} + {Q^} 


( 4 ) 


where 


[C] a S 5 rmmetric matrix of heat capacitance 
[K] a symmetric matrix of heat conduction 

{Q^} a vector of applied thermal loads that are allowed to be time-dependent 
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{Q^} a vector of non-lineai’ thermal loads that depend on temperature 
{ T } a vector of temperatures at grid points 
Equation (4) can be abridged as 


lL] {T} = {Q} 


( 5 ) 


where [ L] , a nonlinear differential operator, is obtained by a comparison of equa- 
tion (4) with equation (5), that gives 


[L] 


[C] • ^ + [K] 
d t 


( 6 ) 


and 


iQ} {Qj} + {Q,} 


(7) 


Constraining a point to a constant temperature, we need to create one new un- 
known and one new equation. The new unknown representing the heat of constraint 
required to sustain the temperature to a prescribed level is defined as 


Tjooi “ heat of constraint (8) 

Adding the quantity T 2 ooj into equation (5) yields the following modified form 


^ 0 ^ 
0 


[L] {T} = {Q} W 


‘200 1 
V 0 J 


(9) 


row for the point 200 


The new equation is defined as 


T 


200 


d200 1 


( 10 ) 
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This quantity q 2 ooi is an additional component in the vector {Q} to be arranged in a 
row corresponding to that of the point 2001, The above equations can be merged to 
become 
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The Bulk Data cards used in modeling are given in figure 5, Case (4). An 
EPOIjNT is used to accommodate the new unknown Tjoo i> Iwo added quantities in 
the merged L-matrix are input through the DMIG cards. The prescribed temperature 
is entered by DAREA with its true value. As usual, a TLQAD2 again serves to refer- 
ence the DAREA and provides an identification to be selected by the DLOAD card in 
the Case Control Deck to activate transient thermal loading conditions. 

The use of DMIG for inputing components directly to the conduction matrix 
requires its function to be selected in the Case Control Deck by a K2PP card. There- 
fore, K2PP= ri-l XEl) , was added. 

Since the modified L-matrix appearing on the left-hand side of equation (11) is 
asymmetric, an altering of the program code to the APP HEAT, SOL 9 is required. 
The DMAP- ALTER packet which was inserted in the Executive Control Deck is repro- 
duced below for reference. 


$ DMAP-ALTER FOR APP HEAT SOL 9 TO FACILITATE ASYMMETRIC CONDUCTION MATRIX 
ALTER 53,53 $ FOR APP HEAT, SOL 9 TO ACCOMMODATE ASYMMETRIC MATRIX HKGG 
MCE2 HUSET,HGM,HRGG,HBGG, ,/HRNN,HBNN,, $ 

PARAM //C,N,PREC/V,N,PREC $ 

UPARTN HUSET,HKGG/HKNNB,HKMN,HKNM,HKMM/C,N,G/C,N,N/C,N,M $ 

MPYAD HKMM,HGM,HKMN /HKM /C , N , 0/C , N , 1/ C , N , 1 / V , N , PREC $ 

MPYAD HKNM,HGM,HKNNB/HKN /C ,N , 0/C , N , 1 /C , N , 1 / V, N , PREC $ 

MPYAD HGM ,HKM,HKN /HKNN/C , N , 1 /C , N , 1/ C , N , 1 / V ,N, PREC $ 

ENOALTER 
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Results of this method are identical to those of Case (1). In addition, the "Tem- 
perature Vector" output for the point 2001 should be interpreted as the heat of constraint 
(in the imits of power, w) as defined by equation (10), 

While the advantages of this modeling technique are self-evident, the disadvan- 
tages are: (A) An increase of the problem size, and (B) The nonstandard definition of 
output. 


For a Time- Varying Temperature Function 

(5) To specify a time-varying temperature function via CELAS2 . QVECT , TABT.FPl , 

TLOAD2 , etc. 

A prescribed temperature at a point varied as a time function can be modeled by 
an extension of the basic technique shown in Case (1). The technique used is to replace 
DAREA with other admissible thermal load cards which would provide with the time- 
varying thermal power capability. QVECT permits its vector components of the direc- 
tional thermal flux to be referenced by the TABLEDl cards where tabulated values of 
temperature vs. time fimctions can be entered. 

The CELAS2 and TLQAD2 serve the same functions as described in Case (1). 

Since QVECT can only input a thermal flux to a point through a boundary surface 
element, CHBDY is therefore required. To define the surface area of the provided 
CHBDY element as well as to specify the surface absorptivity to absorb the radiant 
flux, PHBDY must be included. Moreover, CHBDY can only be referenced by a GRID , 
which, therefore, must replace the SPOINT 200 in the packet of the Bulk Data cards as 
shown in figure 6, Case (5). 

(6) To specify a time-varying temperature function via CELAS2 , QVECT , TLOADl , 

TABLEDl , etc. 

A variation of the modeling technique used in the last example is demonstrated. 
Although TLOAD2 is the only card that must be replaced, in this case by TLOADl , the 
sequence of referencing from card to card to perform specific fmictions has altered 
substantially. The time-varying temperature function does not rely on the provision of 
the vector components of QVECT. Instead, QVECT serves only to facilitate a large 
fictitious thermal load, while tabulated values of a time-varying temperature function 
are entered in TABLEDl which, in turn, is referenced by TLOADl. The packet of Bulk 
Data cards in figure 6, Case (6) shows the detail of inter-relationships among individual 
cards. 
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(7) To specify a time-varying temperature function via CELAS2, QHBDY, TLOADl, 

TABLEDl . etc. 

The preceding case has shown that a temperature vs. time function can be modeled 
via TLOADl and TABLED1» and that a thermal load card is used merely to accommo- 
data a constant value of thermal flux. In the present case, QHBDY is substituted for 
QVECT . Nothing else is altered and the packet of Bulk Data cards is shown in figure 6, 
Case (7). 

(8) To specify a time-varying temperature function via EPOINT, DMIG, DAREA, 

TLOADl and TABLEDl 

As evidenced in Cases (6) and (7), a combination of TLOADl and TABT-EDI is 
well suited to accommodate temperature vs. time functions. 

Therefore, the modeling techniques previously devised for a constant prescribed 
temperature can be modified easily if TLOAD2 , which was used originally in the con- 
stant prescribed temperature case, is replaced by a combination of TLOADl and 
TABLEDl . An application of this substitution to Case (4) has successfully demonstrated 
another method available to the users. The vmiqueness and features of this modeling 
technique are identical to that of Case (4). A packet of relevant Bulk Data cards is 
shown in figure 6, Case (8). 

Details of the illustrated cases with emphasis on the inter-relationships of the 
cards in the Bulk Data Deck and those between the Case Control and the Bulk Data 
Decks together with results are docvimented separately (ref. 5). 


EFFECT OF THE NODE VALENCY 


A node valency is the pattern of element arrangement in an assembled model in 
relation to an associated grid point. Since a number of factors could influence the 
solution accuracy of a finite- element model (ref. 6), the effect of the node valency can 
be controlled by modeling and is, therefore, considered and discussed herewith. 

The input energy to an element, whether from an external thermal flux to the 
surface of an element or from an internal volumetric heating source, is evenly dis- 
tributed to the vertices of the element according to the solution algorithm currently 
implemented in the NTA program. The effect of a node valency on the total energy dis- 
tribution to grid points is essential to the eventual temperature results. As an illustra- 
tion, witiiout computation, six types of node valency available to model a constant area 
A of a rectangle are shown in figure 7. Each rectangular surface is assumed to be cut 
out from an assembled finite-element model along a boundary. The sxirface is heated 
by a constant uniform heat flux q. Then Q. is the amoimt of energy being distributed to 
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that grid point i from all connected elements that define the node valency at i. For 
example, the rectangular surface is modeled by four equal quadrilateral elements as 
shown in figure 7(a). The distribution of energies at three noted grid points are Q j = 
(l/4)Aq, Q 2 =(l/l 6 )Aq and Q 3 =(l/ 8 )Aq. 

When the same rectangular surface is modeled differently by using triangular 
elements as shown in Cases (b) through (f) in figure 7, (i = 1, 2 , 3, 4, wherever ap- 
plicable) of each case shows that at the same location of the centers of these rectangles, 
Qx has varied from 1/3 to 1/6, Q 2 from 1/6 to 1/24, and Q 3 from 1/2 to 1/8 imits. 

This is an indication of imeven distribution of a uniformly applied heat flux to a surface 
of the same area attributed to the node valency. In practice, a structure represented by 
a finite- element model usually employs sufficiently large number of elements arranged 
in a very regular pattern or patterns. This fact implies that those grid points, being 
fully surroxmded by connected elements, would receive the same fraction of energy as 
long as the grid points have the node valency which would result in identical values of 
Qj. It is also seen that the grid points that are isolated or located along boundary lines, 
where the grid points are not fully surrounded by adjacent elements, would be affected 
by the node valency more than their counterparts located on the interior of the surface. 
Therefore, to delineate temperature solution inaccuracy, when the existing software of 
the NTA is used, the selection of uniformly sized and the similar type of elements 
arranged in a regular pattern for interior points is recommended. Using finer discre- 
tization or smaller elements along the boundary lines is also desirable but the engi- 
neering practice must trade-off with other considerations. 

As a last resort, application of thermal loads manually to those isolated grid 
points or those along the boundary lines is a practical method of compensating for 
uneven energy distribution. The Bulk Data card QHBDY is employable for such an 
application. 
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TABLE 1. - THERMO-PHYSICAL PROPERTIES AND DIMENSIONS 


Layer 

Metallic Slab (a) 

Insulation (b) 

Thermal conductivity 
k (W/cm-°C) 

2.25 X 10-3 

0.2 X 10-3 

Thermal capacitance 
PCp (Cal/cm^-°C) 

1.44 X 10-1 

1.163 

Thickness 

A X (cm) 

0.30 

0.40 
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STRUCTURAL SYSTEM 


THERMAL SYSTEM 





Figure 2. The composite slab and its finite-element representation. 
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Figure 3. Time-dependent prescribed temperatures of the hot flowing gas 
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$ BULK DATA DECK 

$ GRID POINTS 

GRID 1 

GRID 2 

GRID 3 

GRID 4 

GRID 5 

GRID 6 

GRID 7 

GRID 8 

SPOINT 200 

$ CONDUCTION ELEMENTS 


. 4 

.. 5 

• • 

0.00 

0.0 

0.0 

0.05 

0.0 

0.0 

0.10 

0.0 

0.0 

0.15 

0.0 

0.0 

0.20 

0.0 

0.0 

0.25 

0.0 

0.0 

0.30 

0.0 

0.0 

0.70 

0.0 

0.0 

1 

2 

2 

3 

4 

4 

5 

6 

6 


0.0022S 0.144 


1101 
1102 
0.000 2 
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1.0 

1.163 


CONVECTIVE BOUNDARY IN CONTACT WITH HOT GAS 


CHBDY 1001 

ei 200 

PHBDY 1002 
MAT4 1003 


1003 
0.07 1 


RADIATION TO AN ENVIRONMENT OE ZERO DEGREE 
DY 9001 9002 POINT 8 


CHBDY 

6CY001 

9001 

9002 1 

PHBDY 

9002 


RADLST 

9001 


RADMTX 

1 

1.0 

PARAM 

SIGMA 

5.676-12 

PAR AM 

TABS 

273. 16 

S INITIAL COMOITIONS 

TEMP 

701 

200 


TEMPO 701 20.0 

$ TIME INCREMENTS AND INTEGRATION STEPS 

TSTEP 601 15 1.0 2 ^ 

63 15 2.00 2 

$ ESTIMATED TEMPERATURES 
TEMPO 1 1927.0 


ENDOATA 


Figure 4. A listing of input data deck of the NTA model. 
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LABEL= FIXED TEMPERATURE AT SPUINT 200. VIA CELA S2 » DARE A f; TL0A02 
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Figure 5. Bulk Data cards used in various modeling techniques for specifying constant prescribed 
temperatures. 
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Figure 6. Bulk Data cards used in various modeling tedhniques for specifying 
time-varying prescribed temperature functions. 
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Figure 7. Typical patterns of the node valency of a constant Area A. 



